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Guidance  law  of  a  homing  missile  is  implemented  using  an  extended  high  gain  observer  without  the 
model  assumptions  on  target  maneuvers.  First,  for  a  class  of  multi-input-multi-output  nonlinear  systems, 
extended  high-gain  observers  are  used  for  output  feedback  control  with  partial  practical  stabilization. 
Then  the  method  is  applied  to  the  guidance  law  implementation  of  a  homing  missile  with  bearing  only 
measurement.  Simulation  results  show  satisfactory  performance. 
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1.  Introduction 

In  the  terminal  guidance  phase  of  a  homing  missile  with  bear¬ 
ing  only  measurement  (BOM)  system,  as  in  the  case  of  a  homing 
missile  with  a  passive  infrared  seeker,  only  line-of-sight  (LOS)  an¬ 
gle  signals  are  available  for  guidance.  The  widely  used  Proportional 
Navigation  Guidance  (PNG)  law,  as  well  as  many  other  improved 
guidance  laws  based  on  PNG,  feeds  back  the  LOS  rate  signals  to 
stabilize  the  LOS  rate.  Thus  to  implement  the  guidance  laws,  it  is 
necessary  to  estimate  the  differentials  of  LOS  angle  signals.  Beside 
the  LOS  angle  and  LOS  rate  signals,  modern  guidance  laws  utilize 
the  relative  range,  relative  range  rate,  and  target  maneuver  signals 
to  cancel  out  the  nonlinearities  and  compensate  for  the  target  ma¬ 
neuvers.  Conventional  implementation  methods  include  observers 
and  digital  filters,  both  of  which  require  the  observability  of  the 
signals  to  be  estimated  [1-3]. 

Difficulties  arise  from  the  lack  of  observability  of  the  signals 
to  be  estimated.  Simple  analysis  shows  that  if  no  additional  mea¬ 
surement  is  available,  different  ranges  and  target  maneuvers  are 
indistinguishable  from  the  LOS  angles  [2],  A  common  approach  to 
tackle  this  difficulty  is  to  incorporate  the  target  motion  dynamics 
to  the  guidance  model  [1].  The  guidance  performance  will  be  seri¬ 
ously  deteriorated  if  there  exists  significant  difference  between  the 
real  target  maneuver  and  model  prediction. 
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In  this  paper,  a  different  idea  is  proposed.  All  the  signals  nec¬ 
essary  to  the  guidance  law  implementation  are  estimated  using  a 
high  gain  observer.  In  the  relative  motion  equations  of  the  mis¬ 
sile  and  its  target,  the  nonlinear  terms,  consisting  of  the  signals 
of  the  LOS  rates,  relative  range  and  its  rate,  are  treated  as  model 
uncertainties,  and  the  target  maneuvers  are  treated  as  external  dis¬ 
turbances.  The  model  uncertainties  and  external  disturbances  can 
be  estimated  by  extending  the  order  of  the  observer.  Using  the 
estimates  of  the  extended  observer,  the  cancellation  of  nonlinear 
terms  and  compensation  for  target  maneuvers  in  the  guidance  laws 
can  be  implemented.  In  doing  so,  the  observability  assumptions  of 
target  maneuver,  relative  range  and  relative  range  rate,  required 
by  other  methods,  are  removed,  and  no  target  maneuver  model  is 
needed  for  the  guidance  law  implementation. 

The  use  of  high-gain  observers  is  a  practical  approach  to  im¬ 
plementation  of  state  feedback  control  of  nonlinear  systems  [4], 
A  high-gain  observer  can  robustly  estimate  the  derivatives  of  the 
output,  which  are  used  to  replace  the  true  states  in  state  feedback 
control.  The  observer  state  will  converge  very  fast,  but  with  a  peak¬ 
ing  phenomenon.  To  protect  the  state  of  the  system  from  peaking, 
the  state  feedback  control  is  required  to  be  globally  bounded.  This 
requirement  is  usually  met  by  saturating  a  continuous  state  feed¬ 
back  function  outside  a  compact  region  of  interest.  The  theoretical 
base  for  the  use  of  high-gain  observers  in  output  feedback  control 
is  the  so-called  separation  principle.  For  the  closed-loop  system 
of  a  class  of  nonlinear  systems  with  output  feedback  control  us¬ 
ing  high-gain  observers,  a  fairly  general  separation  principle  can 
be  guaranteed  in  the  sense  that  the  performance  of  a  globally 
bounded  partial  state  feedback  control  can  be  recovered,  and  the 
performance  recovery  includes  asymptotic  stability  of  the  origin, 
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the  region  of  attraction  and  trajectories  [5],  Theoretical  investi¬ 
gations  and  practical  applications  of  high-gain  observers  include 
stabilization,  nonlinear  servomechanisms,  adaptive  control,  sliding 
mode  control,  robustness  to  fast  un-modeled  dynamics,  discrete¬ 
time  implementation  and  application  to  induction  motors  [6], 

Although  widely  discussed,  literature  concerning  high-gain  ob¬ 
server  technique  includes  only  the  systems  whose  stability  is  with 
respect  to  all  state  variables.  In  [5],  asymptotic  stability  of  the  ori¬ 
gin  is  discussed.  In  [7],  bounded-input-bounded-state  stability  is 
assumed  on  the  zero  dynamics.  In  other  cases,  similar  stability  as¬ 
sumptions,  such  as  input-to-state  stability,  on  zero  dynamics  are 
required.  However,  for  the  guidance  problem  of  a  homing  missile 
against  a  target,  only  the  LOS  rate  signals,  which  are  only  a  part  of 
the  whole  system’s  state  variables,  are  to  be  stabilized  [8-11],  In 
many  engineering  applications,  consideration  of  stability  with  re¬ 
spect  to  part  of  the  system’s  state  variables,  the  so-called  partial 
stability,  is  also  necessary  [12,13],  Other  examples  of  partial  stabil¬ 
ity  can  be  found  in  many  fields  such  as  electro-magnetics,  inertial 
navigation  systems,  spacecraft  stabilization  via  flywheels,  biocenol- 
ogy  [14]. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2,  the 
guidance  dynamics  of  a  homing  missile  are  discussed.  Section  3 
designs  an  extended  observer  to  recover  the  performance  of  a  state 
feedback  control  law.  In  Section  4,  the  extended  high-gain  observer 
is  used  to  implement  the  guidance  law  of  a  missile  in  a  terminal 
guidance  phase,  and  simulation  is  conducted.  We  concludes  the 
paper  with  Section  5. 

2.  Guidance  dynamics  of  a  homing  missile 

For  a  homing  missile  in  the  terminal  guidance  phase,  the  task  of 
guidance  laws  is  to  steer  the  missile  to  intercept  its  designated  tar¬ 
get  with  a  tolerable  miss  distance.  The  general  three-dimensional 
attack  geometry,  in  an  inertial  coordinates  frame,  is  shown  in 
Fig.  1.  M  and  T  represent  the  missile  and  the  target,  respectively, 
both  of  which  are  considered  as  mass  points  in  designing  guidance 
laws,  r  is  the  relative  range  between  the  missile  and  the  target, 
and  vm  and  vt  are  velocity  vectors  of  the  missile  and  the  target, 
respectively.  The  line  connecting  M  and  T  is  LOS,  the  orientation 
of  which  in  coordinates  frame  Ogxgygzg  determines  qp  and  q£, 
the  LOS  angles  in  horizontal  and  longitudinal  planes,  respectively. 

Three-dimensional  relative  motion  between  the  missile  and  its 
target,  as  shown  in  Fig.  1,  can  be  described  as  follows  [18]: 

f  =  rqE  +  rq2^  cos2  qE  (1) 


h  =  ~  jQP  +  2qEqp  tanqg  +  (2) 

qs  = - qE  -qpsinqgcosqp -  (3) 

where  aty  and  amy  are  the  acceleration  components  of  the  target 
and  the  missile  in  the  longitudinal  plane,  respectively,  both  per¬ 
pendicular  to  the  LOS.  dtz  and  amz  are  the  acceleration  components 
of  the  target  and  the  missile  in  the  horizontal  plane,  respectively, 
both  perpendicular  to  the  LOS. 

Guidance  laws  are  designed  to  issue  lateral  acceleration  com¬ 
mands,  amyc  and  dmzc.  The  autopilot  of  the  missile  receives  lat¬ 
eral  acceleration  commands  from  the  guidance  system  and  causes 
aerodynamic  surfaces  to  move  so  as  to  attain  these  commanded 
accelerations.  In  a  guidance  law  design  routine,  the  dynamics  of 
the  missile  autopilot  control  loop  are  often  omitted  to  simplify  the 
design.  Thus,  we  assume 

Qmy  =  Omyc  >  amz  =  Qmzc 

The  aerodynamic  acceleration  capability  of  a  missile  is  physically 
limited  by  such  factors  as  the  existence  of  stall  angle  of  attack, 
the  maximum  defection  of  movable  control  surfaces.  In  an  actual 
missile,  the  airframe  acceleration  must  also  be  limited  in  order  to 
prevent  structural  failure.  Thus  we  have  the  following  constraint 

|Qmyl^QM>  l^mzl  ^  Qm  (4) 

where  om  is  the  maximum  aerodynamic  acceleration  the  missile 
can  provide  perpendicular  to  LOS,  in  both  horizontal  and  vertical 
planes. 

A  missile  with  passive  BOM  system  is  considered.  The  seeker, 
which  has  a  minimum  effective  operation  range,  denoted  by  te, 
provides  the  LOS  angle  measurement  for  the  implementation  of 
guidance  law,  i.e.,  only  LOS  angle  signals,  qp  and  qE,  are  available 
for  guidance  laws  implementation.  In  the  missile’s  engagement 
against  the  target,  if  the  relative  range  r  is  equal  to  or  less  than  te, 
effective  guidance  signals  are  no  longer  available,  and  seeker  gives 
a  command  to  stop  the  guidance.  In  the  remaining  time  of  the  en¬ 
gagement,  the  missile  is  in  a  free  flight  without  guidance.  Note 
that  the  existence  of  te  rules  out  the  singularity  in  the  analysis  of 
guidance  loop.  In  the  free  flight  phase,  the  minimum  range  of  the 
missile  with  respect  to  the  target,  rmin,  is  the  miss  distance,  a  main 
performance  measure  to  be  minimized  by  guidance  law  design. 

In  designing  guidance  laws,  the  LOS  rates  in  both  horizontal 
and  longitudinal  planes,  qp  and  qe,  are  to  be  stabilized.  In  an  ideal 
interception,  i.e.  when  qp=0  and  q£  —  0,  the  relative  motion  equa¬ 
tions  (1),  (2),  (3)  reduce  to 

r  =  0 

which  means  that  the  missile  flies  directly  to  the  target  with  a 
constant  relative  velocity  r  <  0,  and  a  zero  miss  distance  is  to  be 
achieved.  In  an  actual  engagement,  the  miss  distance,  rmjn,  can  be 
achieved  to  a  small  required  value  if  the  LOS  rates  are  kept  small. 
PNG  laws,  a  widely  used  class  of  guidance  laws,  use  feedback  con¬ 
trol,  proportional  to  the  LOS  rates,  to  stabilize  the  LOS  rates.  It  is 
shown  that  a  PNG  law,  with  an  appropriately  selected  navigation 
ratio,  is  optimal  with  respect  to  a  quadratic  cost,  consisting  miss 
distance  and  lateral  accelerations,  against  a  target  without  maneu¬ 
ver  escapes  [10],  The  guidance  performance  deteriorates  with  the 
presence  of  target  maneuvers  [15].  To  keep  the  optimality  in  a  sce¬ 
nario  with  target  maneuvers,  the  PNG  laws  must  be  augmented  to 
include  target  maneuver  compensations.  An  augmented  PNG  law 
takes  the  following  form 

amy  =  Nyqpr  cos qE  +  aty  -  2 rqp  cosqE  +  2rqEqp  sinq£  (5) 
amz  =  -NzqEr  +  atz  +  2  rqE  +  rqj  sin  qE  cos  qp  (6) 
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where  Ny  >  0  and  Nz>  0  are  constants.  The  relative  motion  and 
the  augmented  PNG  constitute  the  closed-loop  guidance  loop.  Note 
that  only  partial  state  variables,  i.e.  qp  and  qE,  are  to  be  stabilized 
in  the  guidance  loop,  and  only  qp  and  qE  are  available  from  mea¬ 
surement  for  feedback  in  (5)  and  (6).  The  other  signals  in  (5)  and 

(6),  unavailable  from  measurement,  are  to  be  estimated  using  an 
extended  high-gain  observer  described  in  the  next  section. 

3.  High-gain  observer  design  and  performance  recovery 

The  three-dimensional  relative  motion  dynamics  developed  in 
the  above  section  can  be  formulated  as  a  general  multi-input- 
multi-output  nonlinear  system  in  the  following  normal  form 


shown  that  for  any  compact  sets  Sz  C  Rm,  S°  C  Rm,  with  Sz  in 
the  interior  of  Sz,  and  Sx  C  Rp  with  Sz  x  Sx  in  the  domain  of 
/o(.,-),  fo  is  continuous  on  Sz  x  Sx,  and  for  any  z( 0)  e  Sz  and 
x(t)  e  Sx,  there  exists  T  >  0,  such  that  z(t)  e  Sz,  Vt  e  [0,  T],  This 
property  is  less  restrictive  than  various  stability  assumptions  [7,5] 
which  can  be  viewed  as  special  cases  of  the  property  where  T  = 
oo.  In  addition,  when  T  is  finite,  unstable  cases  which  admit  no 
finite  escape  times  can  also  be  covered  by  the  property. 

The  goal  is  to  design  an  output  feedback  controller  to  regulate 
the  partial  state,  w,  to  zero,  for  all  initial  states  in  a  given  compact 
set,  while  meeting  certain  requirements  on  the  transient  response. 
For  system  (7)— (10),  if  only  the  selected  part  of  the  state,  w,  is 
considered,  then  for  the  dynamics  of  w,  we  have 


x  =  Ax  +  B[b(x,  z,  d)  +  a(x,  z,  d)u] 
z  =  /o(x,z) 
y  =  Cx 


(7)  w  =  EAx+ EB[b(x,z,d) +a(x,z,d)u] 

(8)  =  Aw  +  B[b(x,  z,  d)  +  a(x,  z,  d)u ]  (11) 

(9)  where  the  n  x  n  matrix  A  and  n  x  p  matrix  B  are  given  by 


where  x  e  Rp  and  zeR”  are  the  state  variables,  u  e  Rp  is  the 
control  input,  d  e  Rp  is  a  vector  of  unknown  external  disturbances, 
and  y  e  Rp  is  the  measured  output.  The  p  x  p  matrix  A,  the  p  x  p 
matrix  B,  and  the  p  x  p  matrix  C,  given  by 


A  =  block diagfAi , . . . ,  Ap], 


'0  1  0  •••  O' 

0  0  1  -0 

0  0  •••  0  1 

_0  0  0  0_ 

'O' 

0 


B  =  block  diagfBi, . . . ,  Bp], 


C  =  block diagfCi, . . . ,  Cp], 


Q  =  [  1  0  •••  0  0 


Ixpi 


where  1  ^  i  <  p  and  p  =  p\- \ - 1-  pp  represent  p  chains  of  inte¬ 

grators.  b(-,  ■ ,  )  and  a(-,  ■  ,■)  are  p-dimensional  nonlinear  function 
vector  and  pxp  nonlinear  function  matrix,  respectively.  The  vec¬ 
tor  weR”  is  a  part  of  the  state  extracted  from  x  by  the  selection 
matrix 


E  =  block  diagfEi , . . . ,  Ep] 

Ei  =  [ ePj-„i+i  ■■■  ePt  ] T  e R"1' xpi ,  1  < i  <  p 

where  n  =  ni  4 - h  np  and  e,  is  a  vector  with  all  elements  be¬ 

ing  zeros  except  the  ith  element  being  1.  After  the  extraction,  the 
remained  partial  state  of  x  is  denoted  by  w.  It  is  assumed  that 


A  =  E  A  =  block  diag[ A  \ , . 
'0  1  0  •••  0 

00  1  •••  0 

Ai=  :  :  •.  : 

0  0  •••  0  1 

_0  0  •••  0  0 


^p] 


B  =  EB  =  block  diagfBi,  ■  ■  • ,  Bp], 


'O' 

0 


Kp 


Since  Eq.  (11)  is  dependent  on  w  and  z,  Eq.  (8)  and  the  dynam¬ 
ics  of  w  cannot  be  dropped  out  of  consideration.  Obviously,  had 
(x,  z,  d)  been  available  for  feedback  and  the  functions  a(-,  ■ ,  )  and 
b(-,  ■  ,■)  been  known,  we  could  have  chosen  a  static  control  law  as 
the  following  feedback  linearization 


u  =a-1(x,  z,  d)(-b(x,  z,  d)  +  v)  (12) 

and  the  dynamics  of  w  would  have  been  reduced  to 


w  =  Aw  +  Bv  (13) 

The  auxiliary  control  v  could  then  be  designed,  for  example,  as  the 
linear  partial  state  feedback  control  v  =  —Kw,  where  (A  -  BI<)  is 
Hurwitz  and  the  solution  of 


w  =  (A  —  BK)w  (14) 

meets  certain  transient  response  requirements.  The  corresponding 
x  dynamics  are 


Assumption  1.  Both  d(t)  and  d(t)  are  bounded  and  d(t )  e  W  c  Rq, 
where  W  is  a  compact  set. 

Assumption  2.  /o  (■ ,  •)  is  locally  Lipschitz  and  continuous  on  its  do¬ 
main,  and  the  nonlinear  function  vector  b(-,  ■ ,  )  and  matrix  a(-,  ■  ,•) 
are  continuously  differentiable  with  locally  Lipschitz  derivatives. 
The  zero  dynamics  equation  (8)  has  no  finite  escape  time. 

Remark  1.  For  the  system  (7)-(10),  suppose  that  only  w  needs  to 
be  stabilized,  i.e.  only  the  stability  with  respect  to  the  set  {w  = 
0}  is  required.  For  the  partial  state  z,  by  Assumption  2  it  can  be 


x=(A-BKE)x  (15) 

By  choosing  the  feedback  gain  matrix  K,  the  partial  state  w  can 
be  regulated  arbitrarily  small,  i.e.  for  any  8  >  0,  there  exists  certain 
T*  >  0,  such  that  ||w(t)||  <  8/2  for  T*  <  t  <  T.  We  call  the  feed¬ 
back  gain  matrix  K  chosen  in  this  way  acceptable.  In  this  case,  for 
the  partial  stability  of  w,  we  could  choose  Vm(w)  —  wTPww,  in¬ 
dependent  of  w  and  z,  where  Pw  =  P*  >  0  is  the  solution  of  the 
Lyapunov  equation 

PW(A  -  BK )  +  (A  —  BK)tPw  =  — Qw 
for  some  Q,w  =  Q. p  >  0. 
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Remark  2.  For  any  c  >  0,  the  set  Qc  =  {Vw(w)  ^  c}  x  R'0-"  x  Rm  is 
positively  invariant  for  the  closed-loop  system.  For  a  given  compact 
set  xS°- xSz°c  S2C,  if  (w(0),  w(0),  z(0))  e  S°  x  S%  x  S°,  then 
by  the  dependence  of  w  on  w  (see  Eq.  (7))  and  Remark  1,  for  any 
T  <  T,  there  are  compact  sets  SL  and  S'z,  dependent  on  c,  V  and 
S°wxS%x.  S°,  such  that  (w,  w,  z)eQ'c  =  (Vw(w)  <  c}  x  S„  x  Sz 
for  0  ^  t  <  V.  The  existence  of  S'w  is  guaranteed  by  the  fact  that 
the  solutions  of  (15)  have  no  finite  escape  time. 

Only  the  measured  output  y  is  available  for  feedback.  Since 
(A,  C)  is  observable  by  definition,  x  can  be  estimated  using  an 
observer.  But  if  linearization  feedback  (12)  is  applied  to  system 
(7)-(10),  then  the  partial  state  z  is  unobservable  from  y.  So  it  is 
impossible  to  estimate  partial  state  z  using  an  observer.  To  imple¬ 
ment  the  feedback  (12),  we  use  the  observer  to  estimate  x,  and  use 
the  nominal  models  of  a(-,  ■ ,  ),  b(-,  ■  ,•)  in  the  control  law.  then  ex¬ 
tend  the  observer  to  compensate  for  the  errors  between  nominal 
models  and  real  functions,  including  the  external  disturbance  d,  as 
discussed  in  next  section. 

The  extended  high-gain  observer,  with  an  additional  integrator 
to  each  of  the  p  chains  of  integrators  in  (7)-(10),  has  the  form 

x  =  Ax  +  B[a  +  b(x)  +  2(x)u]  +  H(e)(y  -  Cx)  (16) 

&=F(e)(y-Cx )  (17) 

where  2(0,  nonsingular,  and  £(■),  are  the  nominal  models  of  o(-, 
■  ,0  and  b(-,  ■ , ■),  and  are  twice  continuously  differentiable  func¬ 
tions,  respectively.  The  extended  observer  state  a  =  [&i, ...  ,&P]T . 
The  observer  gain  matrices  H(e)  and  F(e)  are  given  by 


H(e)  =  block  diag[ffi (e), ....  Hp(e)] 


F(e)  =  diag 


■<±i 

epi+i  ’ 1 


i  =  l . P 


<±il 

ePp+i  J 


where  a\, . . . ,  alp.,a'p.+v  i  —  1, _ p,  are  real  numbers  and  such 

chosen  that  all  the  polynomials 


sA+1  +  c^sA  -| - h a‘p.s  +  ad.+1 ,  i  =  1 , . . . , p 

are  Hurwitz,  and  e  >  0  is  a  small  design  parameter.  For  the  func¬ 
tion  matrix  2(0,  it  is  assumed. 


Assumption  3.  For  all  (w,  w,  z)  e  Q’c  and  d  e  W,  there  exists  a 
constant  Ka,  0  <  Ka  <  1,  such  that 

|  (a(x,  z,  d)  -  a(x))2_1  (x)  ||  <  Ka 

Remark  3.  By  continuity,  all  the  elements  of  o(x,  z,  d)  are  bounded 
over  the  compact  set  12'x  IV.  For  p  =  1,  the  assumption  can  al¬ 
ways  be  satisfied,  as  shown  in  [7],  For  p  >  1,  in  the  special  case 
where  a(x,  z,  d)  is  diagonal,  if  all  the  upper  bounds  of  \aa(x,z,d)\, 
i  —  1, . . . ,  p,  are  known,  then  Assumption  3  can  also  be  satisfied  by 
choosing 


Similar  conclusion  can  be  drawn  in  the  case  where  a(x,  z,  d)  is  di¬ 
agonally  predominant  for  all  (x,  z,  d)  e  S2'c  x  W.  In  other  cases,  it 
must  be  checked  if  there  are  nominal  models  or  not  such  that  the 
assumption  is  satisfied  for  the  given  systems. 

It  can  be  anticipated  that  the  term  H(e)(y  —  Cx)  will  be  0(e) 
after  a  short  transient  period,  thus  Eq.  (16)  can  be  regarded  as  a 
perturbation  of 

x  =  Ax  +  B[a  +  b(x)  +  a(x)u] 

Multiplying  both  sides  by  E,  we  have 

w  =  Aw  +  B[ct +b(x) +  a(x)u]  (18) 

Then  the  following  controller 

u  =a_1(x)(-<7  -  b(x)  -  Kw)  =  i/r(x,  &)  (19) 

where  w  =  Ex,  will  make  (18)  coincide  with  (14).  To  protect  the 
system  from  peaking  in  the  observer’s  transient  response,  as  con¬ 
ventionally  done  in  output  feedback  control  using  high-gain  ob¬ 
servers,  the  control  should  be  saturated  outside  the  region  of  inter¬ 
est,  i.e.  the  compact  set  in  which  the  state  variables  are  contained 
in  the  interval  [0,  T],  For  notation  convenience,  denote  2-1(x)  and 
a_1(x,  z,  d)  by  2(x)  and  a(x,  z,  d),  respectively,  and  let 

I  p 

Mi  >  max  y^a.-Ax, z, d)(— b,(x, z, d)  +  V.) 

(x,z)eS*xSz,deW|-f^  J  x  1  J/| 

1  =  1,  ...,P  (20) 

where  v  =  [vi,...,vp]r  =  —Kw.  Using  the  standard  saturation 
function  sat  [4-6],  we  saturate  the  control  (19)  as  follows 


u  =  Msat(M  V(x, <t)) 

'  Mi  minfl,  |i/q(x,(T)|/Mi}sign(Vn(x,a)) 

_  Mp  min{l,  |  fv(x,  <r)|/Mp}sign(i/rp(x,  ct)) 

where  M  =  diag[Mi , . . . ,  MpJ. 

Let 


(21) 


G(s)  =  diag[Gi(s), 
G;(s)  = 


• ,  Gp(s)] 


“Pi+i 

u\(s)Pi  + 


By  the  diagonal  structure  of  G(s)  we  have  properties  similar  to 
scalar  case  discussed  in  [7]: 


||G(s)  I  =  max  sup  |  Gf(jco)  |  Js  1 

ld<p  CO 

and  when  poles  of  each  Gj(s)  are  assigned  to  be  real  by  appropri¬ 
ately  choosing  the  design  parameters  a’-,  we  have  II G (s) || =  1. 


Denote  the  solutions  of  (14)  and  (15)  by  w*(t)  and  x*(t),  re¬ 
spectively,  and  we  have  the  following  theorem. 


Theorem  1.  For  the  closed-loop  system  (7)-(10),  (16)-(17)  and  (21) 
with  an  acceptable  K  and  Assumptions  1  to  3  satisfied,  if 

||G(es)||ooKa<l  (22) 


sign(2jj(x))  =  sign(cijj(x,  z,  d)) 

|2ii(x)|>  max  {|aii(x,z,d)|},  i  =  l,...,p 

(x,z,d)et?cxW 


then  with  bounded  observer  initial  state  and  each  (w(0),  w(0),  z(0))  6 
sw  x  S^,  x  Sz,  for  any  S  >  0,  there  exist  T*  e  (0,  T)  and  e  >  0  such 
that  for  each  e  e  (0,  e) : 
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||x(t)-x*(t)||  0  ass  -y  0,  uniformly  in  t, 

for  t  e  [0,  T]  with  x(0)  =  x*  (0) 

|  w(t)||  <  S,  for  T*  <t  <T 

The  proof  of  the  theorem  is  given  in  Appendix  A. 


thus  a(-,  ■ ,  )  is  invertible  and  bounded  provided  x3  remains  small. 
By  the  continuity  and  magnitude  limitation  of  target  maneuvers, 
Assumption  1  is  satisfied.  Applying  the  guidance  laws  (5)  and  (6) 
to  the  relative  motion  equations  (l)-(3),  ijp  and  qe  are  stabilized 
and  therefore  x  is  restricted  in  a  compact  set,  denoted  by  Sx  for  1 e 
[0,  T],  For  x  6  Sx  we  can  give  the  description  of  set  Sz  as  follows 


Remark  4.  Better  performance  recovery  can  be  achieved  with  the 
decrease  of  e.  As  pointed  out  in  [5],  the  high-gain  observer,  in 
essence,  estimates  the  derivatives  of  the  state.  Thus  as  shown 
in  [17],  with  the  presence  of  measurement  noise,  the  value  of  e 
suitable  for  practical  use  is  bounded  from  below  due  to  the  noise 
amplification  effect  of  high-gain  observer.  Considering  the  mea¬ 
surement  noise,  the  output  of  the  system  (9)  becomes 

y  =  Cx  +  ns  (23) 

where  ns  e  Rp  is  the  measurement  noise  vector,  bounded  with 
||ns(t)||  <  \i.  In  such  a  case,  the  estimation  error  satisfies  the  ul¬ 
timate  bound 

||x(t)  -  x(t)  |  ^  eci  +  (24) 

for  some  positive  constants  Ci  and  C2.  Since  the  estimation  error 
is  of  the  order  of  magnitude  0(l/er_1),  too  small  s  deteriorates 
the  estimation  performance.  A  tradeoff  must  be  considered  be¬ 
tween  noise  amplification  and  performance  recovery  including  the 
state  reconstruction  speed  and  uncertainties  estimation.  The  in¬ 
equality  (24)  shows  that  the  lower  bound  on  e  is  of  the  order  of 
magnitude  0(/c1/r).  A  switched-gain  observer  scheme  is  proposed 
to  relax  the  tradeoff  in  [17],  to  which  the  interested  reader  is  re¬ 
ferred  for  details. 


4.  Guidance  law  implementation  using  an  extended  high-gain 
observer 


4.1.  High-gain  observer  design  and  guidance  law  implementation 


Take 


'9fi~ 

ip 

Qe 

-Ve  . 

_<Je  _ 

amy 

Qmz 

and  the  relative  motion,  shown  by  Eqs.  (1),  (2)  and  (3),  writes  in 
the  form  of  (7)-(10),  where  p  =  2,  p  =  4,  n  =  2,  m  =  2,  and 


b(x,  z,  d) 


r  bi(x,  z,  d)  1 

1"  -y^*2  +  2x4x2  tanx3  +  1 

[b2(x,  Z,  d)  J 

|_  — ^X4  —  X2  sinx3 cosxi  —  ^  J 

(25) 


a(x,  z,d)  =  [  Z’“SX3  °J 

fo(x,z)  = 

iz 


Z1X3 +ZiX2COSZX3 


(26) 

(27) 


Due  to  the  existence  of  the  effective  operation  range  of  the  seeker 
rE,  in  the  guided  flight  phase  of  the  missile,  we  have  Z\  (t)  ^  rE, 


rm<r^rM,  0  <rm<rE 
rm<r^rM,  rM  <  0 

where  rm,  tm,  rm,  rM  are  parameters  to  describe  all  the  possible  en¬ 
gagement  scenarios  of  the  terminal  guidance.  Denote  by  T  the  first 
time  when  z  goes  out  of  Sz,  then  for  any  t  e  [0,  T],  z(t)  6  Sz.  Since 
r  has  little  change  during  the  terminal  guidance  phase  and  rE>rm, 
so  we  can  appropriately  choose  rm,rM  such  that  there  exists  a  fi¬ 
nite  T  <  T  given  by  r(T)  =  rE. 

In  the  state  space  description  of  the  relative  motion,  target 
maneuver  components,  di  =  aty  and  d2  =  atz,  are  thought  of  as 
unknown  external  disturbances,  and  terms  consisting  of  state  z  are 
thought  of  as  nominal  models  with  uncertain  model  perturbations, 
and  both  uncertainties  are  included  in  the  nonlinear  terms  b(-,  ■  ,■) 
and  a(-,  •  ,-)•  For  this  nonlinear  system,  we  can  design  an  extended 
high-gain  observer  as  follows 

x  =  Ax  +  B[cr  +  b(x)  +  a(x)u]  +  H(e)(y  -  Cx)  (28) 

h  =  F(s)(y  -  Cx)  (29) 

where  x  e  R4  is  the  observer  state  vector,  nominal  models  of 
b(-,  ■  ,•)  and  a(-,  •  ,■), 


are  twice  continuously  differentiable,  a  =  [tf\  <T2]t  is  the  extended 
observer  state  vector,  used  to  estimate  the  uncertainties  including 
both  model  perturbations  and  disturbances,  and  the  observer  gain 
matrices  take  the  form 


T-  0 


Real  numbers  aj,  i  =  1,2,3,  j  =  1,2,  are  chosen  such  that  both 
polynomials 

s3  +  a{s2  +  a]2s  +  aJ3,  j  =  1,2  (30) 

are  Hurwitz,  and  different  values  of  £  >  0  will  be  chosen  in  the 
simulation  to  compare  the  effects  of  state  observation  and  perfor¬ 
mance  recovery. 

Now  we  use  the  estimates  of  the  high-gain  observer  to  imple¬ 
ment  the  guidance  laws  (5)  and  (6).  For  the  nonlinearities  b(-,  ■  ,■) 
and  a(-,  ■ ,  )  in  (25)  and  (26),  state  variable  z3  is  replaced  by  rE,  Z2 
is  replaced  by  a  design  parameter  Vc  <  0  according  to  the  typical 
engagement  scenarios  the  missile  is  designed  to  apply.  As  a  result, 
the  nominal  models  b(-)  and  2(-)  are  as  follows 


a(x)  = 


-y^-x  2  +  2x4X2  tanx3 
—  y^X4  —  x2  sinx3  cosxi 
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The  model  perturbations  b(-,  ■  ,•)  —  &(•).  with  external  disturbances 
d i  and  d2  included,  and  a(-,  ■  ,•)  -  a(-),  are  compensated  for 
with  a.  In  the  guidance  process,  r  <  0  and  T  is  given  implicitly 
by  r(T)  =  rm.  Simple  calculation  can  show  that  Assumption  3  is 
satisfied,  with  Ka  =  1  —  rm/rM  <  1.  When  parameters  at,  i  =  1,2, 
j  ss  1,2,3,  are  such  chosen  that  ||G(s)||oo  =  1,  the  condition  (22) 
is  satisfied.  Thus,  we  have  the  implementation  of  augmented  PNG 
laws  (5)  and  (6)  as  follows 

til  =  Nyx 2Te  cos  £3  —  2VcX2  COSX3  +  2teX4X2  sinx3 

+  r£6-icosx3  (31) 

ti2  =  -NzX4 te  +  2Vcx4  +  tex 2  sinx3 cosxi  -  rE<T2  (32) 

which  are  saturated  as  (21)  with  M\  —  M2  =  cim- 
4.2.  Simulations 

To  compare  the  performance  recovery  of  guidance  law  imple¬ 
mentation  using  the  designed  high-gain  observer,  the  following 
cases  are  studied  in  the  simulations: 

Case  1.  The  augmented  PNG  laws,  (5)  and  (6)  using  state  feedback 
directly  and  assuming  the  target  maneuvers  amy  and  amz  to  be 
known. 

Case  2.  The  augmented  PNG  laws,  (31)  and  (32)  implemented  with 
a  high-gain  observer. 

Case  3.  The  augmented  PNG  laws,  (31)  and  (32)  implemented 
with  a  high-gain  observer  based  on  an  output  with  measurement 
noise  (23). 

Under  PNG  and  augmented  PNG  laws,  homing  missiles  can  con¬ 
duct  interception  in  any  direction  with  respect  to  their  targets  [19, 
20],  So  in  the  simulation  we  choose  an  engagement  scenario  with 
the  following  initial  values  and  target  maneuvers 

r(0)  =  5000  m,  f  (0)  =  -500  m/s,  qp  (0)  =  20° 

qp{  0)  =  0.57°/s,  <?c(0)  =  30°,  qs{  0)  =  2.87°/s 

aty  =40sin^^t+  ^  m/s2,  atz  =  40sin^^t^  m/s2 

We  assume  that  the  missile  is  equipped  with  a  BOM  measurement 
system,  and  the  effective  operation  range  of  the  seeker  is  assumed 
to  be  Te  ~  150  m,  the  maximum  aerodynamic  acceleration  the 
missile  can  provide  is  assumed  to  be  om  =  400  m/s2.  The  design 
parameters  of  augmented  PNG  laws  are  chosen  as  Ny  =  Nz  =  5, 
the  design  parameter  Vc  is  chosen  as  Vc  =  —550  m/s. 

In  the  simulations  for  all  Cases  1  through  3,  we  set  ui  =  0  and 
U2  =  0  when  zi(t)  <  r£.  In  Case  1,  the  miss  distance  is  rmin  = 
1.30  m.  In  Case  2,  we  first  choose  s  =  0.02,  and  the  correspond¬ 
ing  miss  distance  is  rmjn  =  1.49  m,  similar  to  that  in  Case  1.  The 
LOS  rates  X2  =  qp,X4=qe,  and  their  estimates  £2,  £4,  are  shown  in 
Figs.  2  and  3,  respectively. 

From  Figs.  2  and  3,  we  can  see  that  the  LOS  rate  estimates 
converge  rapidly  to  the  real  LOS  rate  signals,  and  simulation  with 
other  values  of  s  also  shows  the  smaller  the  e,  the  faster  the  con¬ 
vergence.  The  results  for  other  values  of  e  are  omitted.  In  the 
simulation  results,  we  can  also  find  the  peaking  of  the  observer 
state,  which  is  a  feature  of  high-gain  observers,  and  the  smaller 
the  e,  the  larger  the  peaking.  But  the  boundedness  of  the  nominal 
models  and  saturation  of  control  protect  the  state  of  the  system 
from  peaking. 


Fig.  2.  LOS  rate  estimation:  qp. 


time(s) 


Fig.  3.  LOS  rate  estimation:  qe. 

To  see  the  performance  recovery  of  the  high-gain  observer,  we 
set  e  =  0.002,  and  compare  amy  and  amz  at  different  values  of  e  in 
Case  2  and  those  in  Case  1,  as  shown  in  Figs.  4  and  5.  It  can  be 
seen  that  with  the  decrease  of  the  value  of  e,  the  acceleration  com¬ 
mands  implemented  with  high-gain  observers  converge  to  those 
using  direct  state  feedback  with  amy  and  amz  known.  This  phe¬ 
nomenon  is  consistent  with  the  dependence  of  the  LOS  rate  con¬ 
vergence  on  the  value  of  e,  as  discussed  above.  When  e  =  0.002, 
the  LOS  rate  signals,  qp,  qE,  are  shown  in  Figs.  6  and  7. 

In  Case  3,  both  components  of  the  measurement  noise  are  as¬ 
sumed  to  be  uniformly  distributed  random  variables  taking  values 
between  —0.6  mrad  and  0.6  mrad,  and  they  change  values  every 
0.0002  s.  When  s  =  0.02,  the  miss  distance  is  rmin  =  1.50  m,  little 
larger  than  that  in  Case  3  with  the  same  value  of  e,  which  means 
that  for  the  high-gain  observer,  the  chosen  e  achieves  a  satisfactory 
tradeoff  between  performance  recovery  and  noise  amplification. 
For  larger  or  smaller  values  of  e,  the  miss  distance  increases  due  to 
the  bad  estimation  or  noise  amplification.  The  real  LOS  rate  signals 
and  their  estimates  are  shown  in  Figs.  8  and  9,  and  the  accelera¬ 
tion  commands  are  shown  in  Fig.  10. 

It  is  worth  noting  that  with  the  presence  of  measurement  noise, 
different  magnitudes  of  measurement  noise  signals  decide  differ¬ 
ent  values  that  e  takes  to  achieve  an  optimal  tradeoff  between 
performance  recovery  and  noise  amplification,  but  the  statistics  of 
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Fig.  5.  Performance  recovery  of  amz. 


Fig.  6.  LOS  rate:  qp. 


Fig.  9.  LOS 
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&e(y) = 


y  for  0  <  y  <  1 

y+(y-l)/£-0.5(y2-l)/£  forl<y<l  +  £ 
l  +  0.5e  fory>l  +  e 


Since  the  function  ge  is  continuously  differentiable,  and 


BKpi+D 

d 

d&  : 

_  rlp(.pp+ 1) . 

=  —I  —  A  “(x,  z,  d) 

x  diag[gg(i/fi(x,  g'e(irp(x,  a )/Mp)]a_1(x) 

It  is  easy  to  show  that  |gg(y)|  <1  for  any  y  e  R,  then  in  view 
of  Assumption  3,  we  can  see  that  the  change  of  variable  (34)  is 
globally  well  defined. 

The  state  equations  of  the  fast  variables  are 


measurement  noise  are  not  necessary  to  the  high-gain  observer 
design,  which  is  different  from  the  filtering  method,  for  which 
imprecise  knowledge  of  the  measurement  noise  statistics  seri¬ 
ously  deteriorates  the  estimation,  even  results  in  instability.  Also 
note  that  in  both  Cases  2  and  3,  no  a  priori  model  assumptions 
are  made  on  target  maneuvers,  although  the  target  maneuvers 
are  unobservable  from  the  BOM  measurement.  While  for  filtering 
method,  either  the  target  maneuver  models  are  assumed  [1],  or 
guidance  laws  are  adapted  to  improve  observability  [2], 


£%  =  -a'jVil  +  WO'+D’  1  <  j  <  Pi  -  1 
erjipi  =  -<xlp.  m  +  mpt+D  +  A°(*> *’  o'. d >  £) 

p 

eriiiPi+V)  =  -otpi+im  -  Aifc(x’ z’  o)/Mk) 

x  aw(x)aj,(+1  rjn  +  £  [  Af  (x,  z,  d) 


5.  Conclusions 

For  a  class  of  multi-input-multi-output  nonlinear  systems,  ex¬ 
tended  high-gain  observers  are  used  in  the  output  feedback  control 
design.  Motivated  by  the  guidance  of  a  homing  missile,  partial 
practical  stabilization  is  considered  for  the  class  of  systems.  Since 
only  a  part  of  state  variables  are  practically  stabilized,  no  non¬ 
minimum  phase  conditions  are  required,  unlike  the  common  cases 
where  stabilization  of  full  state  is  considered.  So  systems  with  un¬ 
stable  zero  dynamics  can  also  be  included  in  this  paper,  only  if 
they  admit  no  finite  escape  times. 

As  done  in  [7],  the  high-gain  observer  is  extended  to  estimate 
the  perturbation  due  to  model  uncertainty  and  disturbance.  The 
performance  recovery  is  proved. 

Then,  the  extended  high-gain  observer  is  applied  to  the  guid¬ 
ance  law  implementation  of  a  homing  missile  with  bearing  only 
measurement.  Simulation  results  show  the  satisfactory  perfor¬ 
mance  recovery  of  the  extended  high-gain  observer. 

Appendix  A.  Proof  of  Theorem  1 

Fast  variables  are  scaled  estimation  errors  defined  as  follows 
=  (33) 

P 

Vi<jH+V)  =  A?(x.  Z,d)  +  J2  A?k(x,  z,  d)Mkge(i/rk(x,  a)/Mk)  -  at 

1  <  i  ^  p  (34) 

where 

A-’fx,  z,  d)  =  bj(x ,  z,  d)  -  bi(x) 

A“k(x,  z,  d)  =  a,k(x,  z,  d)  —  a ;k(x),  1  <  i,  k  <  p 


+  A?k(x,  z,  d)MkgE(^k(x,  a)/Mk ) 

k=l 

+  £  A“k(x,  z,  d)g;  (MX, &)/Mk)  d-^Mlx  1 

k=l  X  J 

for  1  <  i  <  p,  where 
A°(x,  z,  x,  ff,d,  £) 

P 

=  bi(x)  -  bi(x)  +  ^aik(x,  z,  d)Mk[gE(fk(x,  a)/Mk) 

-&(^(x,ff)/Mk)] 

+  Yl[dik(x)Mkgs(irk(x,  a)/Mk) 
k=  1 

-aik(x)Mkge(Mx,t)/Mk)] 

p 

+  ^[a(x,  z,  d)  -  a(x)]Mk[sat(^k(x,  a)/Mk) 
k= 1 

-  ge(fk(X,v)/Mk)] 

The  fast  subsystem  is  given  by 


“pp+i'Jp1 


+  £[Bi  A2  +  B2A3] 


(35) 


and  the  odd  function  ge( •)  is  defined  by  [7] 


where 
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,  Ap],  A\  — 


[[:.] . [:. 

[['•'] . [■•■ 


B2  =  block  diag 
Ai  =  A“(x,  z,  d ) 

x  diag [g'g (irx (x, ct)/Mi) ,  ...,g'E(^p (x, &)/Mp)]d~ 1  (x) 
'  5e(^i(x,  "| 

A2  =  Ab(x,  z,  d)  +  Aa(x,  z,  d)M 


g£(fP(x,a)/M1)_ 


3Vfl(x,o-)  • 


Sg(Vfi(x>5')/Mi) 


A3  =  -A0(x,  Z,  X,  a,d,E) 
s 

By  Assumption  3  and  the  property  that  |g'(y)|  ^  1,  we  can 
show  that  ||A]  ||  s:  1.  By  the  definition  of  variable  change  (33)  and 
the  definition  of  function  ge,  we  can  show  that  A 0/s  is  locally 
Lipschitz.  Assumption  2  ensures  that  A2  and  A3  are  locally  Lip- 
schitz  and,  for  (x,  z)  e  12',  bounded  from  above  by  affine  in  1 1 77 1 1 , 
uniformly  in  e.  Let 

C  =  block diag[a^+1  Ci, . . . ,  appp+,C„\ 

and  without  considering  the  term  e[Bi A2  +  B2A3],  we  first  con¬ 
sider  the  following  system 


sf]  =  A17  —  Bi  AiCty 


a  =  Ah(x,  z,  d)  +  A“(x,  z,  d)M 


ge(irp(x,6)/Mp)  _ 


By  (19),  it  follows  that 

VKx,  a)  =  a_1(x)(-d-  -b(x)  -  Kw) 

From  which  and  (38),  we  can  show  that 
ir(x,a)  =a-1(x, z, d)(-b(x, z, d)  -  Kw) 
holds  up  to  an  0  (e)  error.  Hence,  we  have 
Vr(x,<x)  =  a-1(x, z, d)(-b(x,z,d)  -  Kw)  +  0(e) 

For  (w,  w,  z)  6  Q",  i.e.,  in  the  linear  region  of  saturation  function, 
the  closed-loop  system,  under  feedback  law  u  =  i/r(x,  a),  can  be 
represented  as  follows 


x=(A-BKE)x+0(e) 

z=/o(x,z) 

e?)  =  Ar)  -  BiAiCt/+  0(e) 


(39) 

(40) 

(41) 


After  the  interval  Ti(e),  the  partial  state  rj  will  reach  the  set 
{W (17)  <  e2Ci),  even  if  the  initial  value  17 (Q)  is  outside  the  set. 
Then  for  Ti(e)  <  t  ^  T,  x(t)  satisfies  (39).  If  x*(0)  =  x(0)  then  we 
have 


x(t)  -  x*  (t)  =  0  (Ti  (e)) ,  for  0  <  t  <  Ti  (e) 


(42) 


thus,  x(T  1(e))  —  x*(Ti(e))  =  0(T](e)).  Due  to  the  continuous  de¬ 
pendence  of  the  solutions  of  differential  on  initial  conditions  and 
parameters,  it  is  concluded  that 

x(t)  —  x*  (t)  =  0  (e)  +  0  (Ti  (e)) ,  forTi(£)<t<r  (43) 

Eqs.  (42)  and  (43)  lead  to 


x(t)  -x*(t)  =  0(£)  +  0(Ti(£)),  for  0  <  t  <  T 


(44) 


(36) 


The  system  can  be  viewed  as  a  negative  feedback  connection  of 
a  linear  system  and  nonlinear  gain  Aj,  where  the  transfer  func¬ 
tion  matrix  of  the  linear  system,  C((£s)/  -  zl)-1Bi,  is  C(£s), 
which  is  strictly  proper.  By  Assumption  3,  it  can  be  proved  that 
if  || G (s) 1 00  <  K~\  then  [/  +  KaG(£S)][/  -  KaGffis)]-1  is  strictly 
positive  real,  for  the  details  of  the  proof,  see  Example  7.1  in  [16], 
Then  according  to  the  circle  criterion  of  absolute  stability,  it  can  be 
shown  that  the  system  (36)  will  be  globally  exponentially  stable, 
and  has  a  quadratic  Lyapunov  function  W(t])  —  r)T  P17.  Using  W(rj) 
as  a  Lyapunov  function  candidate  for  (35),  we  can  conclude  that 
for  any  ci  >  0,  17  enters  (W(r/)  ^  e2ci)  in  finite  time  Ti(£),  where 
lim£_>o  Ti  (£)  =  0.  Since  the  control  is  saturated  by  (21),  the  right- 
hand  side  of  (7)  is  bounded  in  [0,  Ti  (£)]  uniformly  in  e,  thus  if  £  is 
small  enough  such  that  Tj(s)  <  T  we  have  (w(t),  w(t),z(t))  e  Q'c 
for  t  6  [0,  Ti(£)]  by  Remark  2.  By  Remark  1,  there  is  a  compact 
set  U”  with  Q’c  in  the  interior  of  it,  such  that  for  V  <  t  <  T, 
(w,  w,  z)  e  i2". 

While  (w,  w,z,  17)  e  12"  x  (W(?j)  <  £2Ci},  we  have  17  =  0(£), 
and  thus 

(37) 


Since  lim^o  Ti(e)  =  0,  we  have  ||x(t)  -x*(t)||  0  as  £  ->  0.  For  a 

given  set  of  initial  conditions  S°  and  a  given  &  >  0,  an  acceptable  K 
can  guarantee  that  ||w*(t)||  =  ||Ex*(t)||  <  8/2,  T*  <t  <  T,  for  some 
time  T*  <  T.  By  (44),  we  can  choose  s  small  enough  such  that 
||x(t)  -  x*(t)||  <  8/2,  thus  ||  w(t)  -  w*(t)||  <  8/2.  Then  we  have 

II  W(t)  |  =  ||w*(t)  +  w(t)- w*(t)|| 

<||w*(t)||  +  ||w(t)-w*(t)|| 

<  ^  +  U  =  8,  for  T*  <  t  <  T 
with  which  the  proof  is  concluded. 
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